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Abstract 

We give a new representation as tempered distribution for the energy-momentum 
tensor of a system of charged point-particles, which is free from divergent self- 
interactions, manifestly Lorentz-invariant and symmetric, and conserved. We present 
a covariant action for this system, that gives rise to the known Lorentz-Dirac equa- 
tions for the particles and entails, via Noether theorem, this energy-momentum 
tensor. Our action is obtained from the standard action for classical Electrody- 
namics, by means of a new Lorentz-invariant regularization procedure, followed 
by a renormalization. The method introduced here extends naturally to charged 
p-branes and arbitrary dimensions. 
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1 Introduction and Summary 

The dynamics of a classical charged point-particle interacting with an electromagnetic 
field is described by the Maxwell and Lorentz equations, 

d^ u = f, = d»A u -d u A^, (1.1) 

^ = eF^(y(s))u u (s). (1.2) 

We parametrize the worldline of the particle y' 1 (s) by the proper time s, and define four- 
velocity and four-acceleration respectively by, 

dy^ „ du^ 
u = -f-, w» = — . 
as as 

The current is, 

^(x) = e y u"(s) 5 4 (x - j/(s)) ds. (1.3) 
Maxwell's equations admit the general solution, 

t - b LW + ^in > 

where is a free external field, and amounts to the retarded Lienard-Wiechert 
(LW) field. 

When substituting the total field strength into the Lorentz equation one faces the 
problem that the LW-field is infinite on the particle's trajectory, and the quantity, 

Flwivis)), (1-4) 

diverges. This infinite self-interaction is, of course, a consequence of the point-like nature 

of the charged particle. 

There are several methods of dealing with this divergence, all of them leading to the 

same conclusions: the quantity (1.4) has a divergent part, which leads to an infinite 

classical renormalization of the particle's mass, and a finite part that amounts to, 

^uv / / \ \ e / ,,dw u „dw^\ . 

^W.»U = - S (^ (1-5) 

Substituting this into the Lorentz-equation one obtains the relativistic Lorentz-Dirac 
equation, 



which has to be considered as the equation of motion for the charged particle, replacing 
the original ill-defined Lorentz-equation (1.2). Here, and in the following, we indicate 
the square of a four- vector with w 2 = w^w^. Strictly speaking this equation can not 
be derived from the fundamental equations (1.1) and (1.2) of classical Electrodynamics 
(ED), but has to be postulated. 

In its relativistic form equation (1.6) has been derived first by Dirac [1] in 1938, in an 
attempt to describe the motion of a classical relativistic charged point-particle, taking 
into account its self-interaction. In this seminal paper the equation was derived using in 
an essential way the principle of energy-momentum conservation. This conservation law 
represents still the basic motivation for the equation itself. 

Despite of the fact that (1.6) is a well-defined equation and ensures energy-momentum 
conservation, it exhibits several unsatisfactory features, both physically and mathemat- 
ically. The key physical drawbacks are that the equation is of the third order in the 
time-derivative, and that it admits a class of solutions, so called runaway solutions - see 
e.g. [2] - which lead to exponentially growing four-velocities, even in absence of external 
forces. When these solutions are eliminated imposing suitable boundary conditions at in- 
finity [3, 4], it turns out that the motions allowed by these conditions show up an acausal 
behavior, in the sense of a preacceleration: the acceleration is felt before the external 
force begins to act. Actually, this preacceleration occurs on a time scale of order, 

e 2 

r = l(T 23 s, 

6nm 

which is of a factor 1/137 shorter than the time scale r q = H/m, where, due to the 
Heisenberg uncertainty principle, quantum effects become relevant and the classical theory 
has to be abandoned. Hence, physically one can interpret this acausal behavior as an 
inconsistency of classical ED of point-particles, which is cured by the quantum formulation 
of the theory, see e.g. [5] 3 . 

Accepting (1.6) as the correct equation of motion for the particle, three questions 

3 One can, actually, maintain a second order Lorentz-equation which is free from singularities, if one 
introduces a finite cutoff of physical origin on the charge distribution of the would-be point-particle. 
This is the philosophy pursued mainly in [6, 7] and references therein. The main unsatisfactory aspects 
of this approach are its lack of universality, and the problems related with the implementation of exact 
energy-momentum conservation. In the present paper we follow Dirac's original proposal and insist on 
point-like charges, and hence on equation (1.6). 
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arise automatically, to which we will give new answers in this paper. The first question 
regards local energy-momentum conservation, in terms of a well-defined conserved energy- 
momentum tensor. The second regards the existence of an action from which (1.6) can 
be derived, and the third question is if the energy-momentum tensor can be derived from 
this action via Noether's theorem. 

Regarding the first question we mention that, although relying on conservation crite- 
rions, in his paper [1] Dirac does not construct an energy-momentum tensor, nor does 
he give a recipe to compute the four-momentum enclosed in a volume V containing the 
particle. The main obstruction to such a construction arises from the particular form of 
the naive energy-momentum tensor of the electromagnetic field, 



where for simplicity we have set the external field to zero. If we indicate with R the 
distance from the particle at a certain instant, the fields F^w exhibit near the particle 
the standard Coulomb-like integrable 1/R 2 singularities, and hence exhibits the non- 
integrable l/R 4 singularities. As a consequence the momentum integrals Py = f v d 3 XT®£ l 
diverge. In a mathematical language the problem can be stated as follows: while the com- 
ponents of F£w are tempered distributions, i.e. elements of S' = S'(R 4 ), the components 
of r^, being products of the formers, are not. 

To obtain a well-defined energy-momentum tensor one has to replace with a 
tensor whose components are elements of S' such that, a) coincides with in 
the complement of the particle's worldline, and b) such that the resulting total energy- 
momentum tensor is conserved, if (1.6) holds. Such a tensor admits then automatically 
finite four-momentum integrals. This program has been accomplished only rather recently 
by Rowe in [8], relying partially on [9]. His energy-momentum tensor is defined, rather 
implicitly, as the distributional derivative of the sum of other elements of S', and entails 
the additional drawback that it is not manifestly symmetric and traceless. On the other 
hand it can be shown that the resulting tensor is uniquely determined by the requirements 
a) and b). 

In this paper we present an alternative and simple representation of this tensor, rely- 
ing on a new Lorentz-invariant regularization procedure, followed by a renormalization, 
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which leads to a manifestly symmetric, traceless and invariant tensor. In particular, our 
procedure follows the physical intuition that the renormalized tensor should differ from 
(1.7) only on the particle's trajectory: correspondingly in our renormalization scheme we 
subtract only (divergent) terms which are supported on the world-line, i.e. terms that 
are proportional to delta-functions supported on the trajectory. 

An alternative approach to four-momentum conservation, based on a different reg- 
ularization procedure, has been proposed in [10]. This approach leads to finite four- 
momentum integrals, but it does not allow to construct an energy-momentum tensor. 
Moreover, it lacks manifest Lorentz invariance. 

The second question opened by (1.6) is if this equation can be derived from an action 
principle. It is immediately seen that the resulting equation can not be derived from a 
canonical action. Indeed, such an action would be necessarily of the type, 



where C se if is an invariant function of the kinematical variables y^, u^ 1 , w^ 1 etc, giving rise 
to the self-interactions at the r.h.s. of (1.6). But for dimensional reasons the unique term 
with the right dimensions would be the invariant C se if oc u^w 1 *, and this is zero. 

Since there exists no canonical action leading to (1.6), a variational principle giving 
rise to this equation bears necessarily some unconventional features. The first proposal 
for an action was made in [11], relying on the previous attempt of [1], and it requires the 
introduction of two gauge fields: a source-free field and a field generated by the particle 
via retarded and advanced potentials. As additional - and more serious - drawback 
this action contains non-local Coulomb interactions between the particles. The more 
recent attempt in [10] provides an action involving a single gauge field, but employs a 
regularization which is not manifestly Lorentz-invariant. 

In this paper we present an action which is obtained from the standard action for 
ED, with a single gauge field, in a very simple and natural manner: we introduce a 
manifestly Lorentz-invariant regularization in the standard action, and renormalize it 
upon subtracting an infinite mass term. The so obtained action gives rise to (1.6) in the 
following way. First one varies the action with respect to the gauge field, obtaining a 
regularized version of the Maxwell equation. Then one varies the action with respect to 
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the coordinates of the particle and substitutes in the resulting equation the solution of 
the regularized Maxwell equation. The limit of this equation, when the regulator goes to 
zero, is (1.6). 

The third question regards the derivation of the energy-momentum tensor through 
Noether's theorem, from a covariant action. As far as we know this problem has been 
addressed only in [12], using a formalism introduced in [13]. This paper is based on a 
manifestly Poincare-invariant, but rather unconventional distribution-valued Lagrangian 
density. Its most unusual characteristics are that the minimal interaction term A^j^ is 
absent, and that the self-interaction of the particle is represented by a total derivative. 
The most unsatisfactory feature of the resulting action is, however, that it does not give 
rise to the equations of motion, i.e. (1.1) and (1.6). Eventually the role of the Lagrangian 
of [12] is restricted to give rise, through translation invariance and Noether's theorem, to 
the conserved energy-momentum tensor given in [8], upon imposing the Lorentz-Dirac 
equation (1.6) "by hand". 

The Lagrangian which we propose in the present paper, being in particular manifestly 
Poincare invariant, allows to apply Noether's theorem in a canonical way, and the resulting 
energy-momentum tensor arises precisely in the new form mentioned above. 

The formalism used in this paper relies on the power of full-fledged distribution theory. 
We remember, indeed, that in the case of point particles the energy-momentum tensor as 
well as the Lagrangian density are necessarily distribution valued. The other key ingre- 
dient is a manifestly Lorentz-invariant regularization procedure. Within this approach, 
in summary, we can set up the entire lagrangian formalism for the ED of point parti- 
cles as specified above: construction of a conserved and well-defined energy-momentum 
tensor, construction of an action which gives to the correct equations of motion, and 
implementation of the Noether-theorem, based on this action. 

The method developed in this paper is based on the invariant regularization of the 
Green function for the D'Alambertian in eq. (2.6), amounting to the replacement x 2 — > 
x 2 — e 2 , that permits a clear and simple distinction between finite and divergent terms in 
the energy-momentum tensor. Since the Green function in arbitrary dimensions depends 
always on x 2 , see [14], our method is very promising for the construction of finite and 
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conserved energy-momentum tensors for charged systems, for which this tensor has never 
been constructed before, [15]. Examples include point particles in a curved background 
[16], systems of dyons [17], point particles in higher dimensions [18], and extended ob- 
jects in arbitrary dimensions [19]. The knowledge of such a tensor allows a systematic 
quantitative analysis of radiation effects in these systems. 

The plan of the paper is the following. In section two we introduce the regularization 
used in the paper. In section three we present and illustrate our main results, i.e. the 
new expression of the renormalized energy-momentum tensor together with its main 
properties, the renormalized Lagrangian density giving rise to the equations of motion, 
and the Noether-theorem. Section four is devoted to the proof of the properties of the 
renormalized energy-momentum tensor, and to a comparison of this tensor with the one 
proposed by Rowe [8]. The variational results are proved in section five. Some technical 
details are relegated to appendices. 

For simplicity we present our results in detail in the case of a single point particle 
and for a vanishing external field . The corresponding results for the general case are 
summarized in section (3.3). 

2 Regularization 

In Lorentz-gauge the Maxwell equations amount to, 



where for a point-particle with charge e and worldline y M (s), the current is given in (1.3). 
The solution of these equations can be obtained in terms of the retarded Green function 
for the D'Alambertian, supported on the forward light cone, 



where x M = (x°,x) are the space-time coordinates, and H denotes the Heaviside step 
function. Omitting the external potential Af n one gets as solution the retarded Lienard- 
Wiechert potential, 



= f, = o, 



(2.1) 





(2.2) 



At 



G*j» = 



e 




(2.3) 



LW — 



4tt (x„ - y u (s))u v (s) 



s=s(x) 
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where * denotes convolution, and the scalar function s(x) is uniquely fixed by the retarded 
time condition, 

( x -y(s)f = 0, x >y°(s), (2.4) 

consequence of the 5-function appearing in (2.2). 

Along the worldline of the particle, i.e. for x^ = y M (A) for some A, the LW-potential 
Al w (y(X)) diverges. In the present notation this follows from s(y(X)) = A. As a conse- 
quence also the LW field-strength, 

diverges along the particle worldline, meaning that a charged particles feels an infinite 
self-interaction. 

To smooth these singularities, while maintaining the structure of the Maxwell equa- 
tions unaltered, we propose to regularize the retarded Green function in a Lorentz- 
invariant way. Choosing an UV-regulator e > 0, with the dimension of length, we 
introduce a regularized Green function [20], 

G e (x) = ±-H(x°)8(x 2 -e 2 ) ee -L 6 + (a*-e*) = \ 5 (x° - ^ + eA , (2.6) 
2tt 2tt Att^/\x\ 2 + e 2 V / 

which is still an element of S'. The support of G £ is given by the positive-time sheet of 
the hyperboloid, 

2 2 

x = e . 

This regularization preserves therefore also causality. 
We can now define a regularized LW potential, 

e u^(s) 



(2-7) 

s=s E (x) 



Air (x v - y u (s))u u (s) 

which differs from (2.3) only through the fact that s(x) is replaced by s £ (x), solution of 
the regularized retarded-time equation, 

(x - y(s)) 2 = e 2 , x°>y°(s). (2.8) 

The potential is a regular field, belonging to C°° = C°°(R 4 ), and in particular it is 
regular on the support of j^, i.e. on the worldline of the particle. Also the regularized 



field strength, 

F£ v = d»A u £ - d u A», (2.9) 

is an element of C°°, and hence everywhere regular. We illustrate these properties giving 
the explicit formulae for a static particle, at rest in x = 0, 

^^^whw' Ae=0, ^ = ^wtw^' 4 = ' (2 ' 10) 

which are all regular in x = 0. 

One of the advantages of this regularization is that it preserves the Lorentz-gauge. 
Indeed, due to the properties of the convolution and thanks to the fact that the current 
is conserved, the definition (2.7) gives directly, 

= 0. 

We can then define also a regularized and conserved current through, 

% = dpFr = nAZ, d u f e =0, (2.11) 

that belongs to C°°, too. 

Since the derivative is a continuous operation in S', it is clear that we have the following 
limits, 

S'-hmA^ = A£ w , (2.12) 
S'-hmFr = Fjjwi (2-13) 
S'-lim^ = f, (2-14) 

where S'— lim denotes the limit in the (weak) topology of S' 4 . Moreover, since by 
construction G £ is Lorentz invariant, F^, as well as are Lorentz covariant tensor 
fields. 

We end this section writing more explicit expressions for F^ v and jg, since these will 
be used in the following. For this purpose we define the vector fields, 

R tM (x) = x t "-y'", ^(x) = {uR)w' 1 -{wR)u' 1 ) (2.15) 



We remember that this means that these limits hold on every test function ip £ S = 5(R 4 ), 

VanA%(<p)=A* w (<p), etc. 
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where the kinematical quantities y, u and w in these expressions are all evaluated at the 
the proper time s e (x), determined from (2.8). In the following this functional dependence 
in the regularized quantities will always be understood. For the scalar product of two 
vectors we use the notation (ab) = a^. In particular we have, see (2.8), 

R^R* = (R ) 2 - \R\ 2 = £ 2 . (2.16) 

The expression for the regularized LW field strength can be derived in the same way as 
the known expression for the standard LW field, (2.5). To calculate it from (2.7) one 
needs the derivative of the function s e (x), that is obtained differentiating (2.16), 

ds £ R„ 



dx» (uR) ' 

The rest of the calculations is a bit lengthy but straightforward. With the above conven- 
tions the results read, 

* = 5R (2 ' 17) 

Fr = ^—^—,{R^u u -R v u^ + R^ u -R u ^), (2.18) 
47T (uR) 3 

« = ^felM^-H.I + ^K + A')), (2.19) 

where we defined, 

V = (2.20) 

as 

The formula (2.18) coincides, actually, with the known expression of the LW field strength 
Flw _ see f° r example [4] - the unique difference being that the kinematical quantities 
are evaluated at the retarded proper time s e (x), rather then at s(x). 

In the following we will use an asymptotic condition on the worldlines y^i^s) of the 
particle. We suppose that the particle's acceleration is zero before a certain instant 5 , 

tu"(s) = 0, for s < s. (2.21) 

Since we have in any case lim^i^oo s e (x°, x) = — oo, this asymptotic condition implies that 

at fixed time the acceleration w^(s £ (x)) in (2.18) vanishes for sufficiently large \x\. The 

5 Actually, the results we obtain in this paper are valid also if the acceleration vanishes sufficiently fast 
for s — > — oo. 



same holds for the corresponding acceleration w^(s(x)) in the unregularized field F{^ L . 
This implies that the regularized and unregularized LW-fields have the same fixed-time 
Coulomb-like asymptotic behavior for large \x\, 

F!r ~w F ™~w (2 - 22) 

These behaviors, that are valid also if the acceleration vanishes sufficiently fast in the 
remote past, will ensure in particular that the total four-momentum of the electromagnetic 
field is finite. 

3 Main results 

3.1 The energy— momentum tensor 

We begin with the construction of a consistent energy-momentum tensor, accepting as 
equation of motion (1.6). 

As explained in the introduction the naive energy-momentum tensor in (1.7) is 
not an element of S 1 - although Fj%y is - due its 1/|-R| 4 behavior near the worldline 6 . 
Therefore its derivatives do not even make sense, and the question what the quantity 
"^ r <fm" amounts to, is meaningless. 

To construct an energy-momentum tensor that is an element of S' we start from the 
regularized LW field (2.18), and define the regularized electromagnetic energy-momentum 
tensor as, 

= F^F ep " + rT F^ a F £pa . (3.1) 

This tensor belongs to C°° and shares with the asymptotic behavior for large \x\ at 
fixed time, implied by (2.22), 

~ — , ~ — . (3.2) 

\x\ 4 \x\ A 



6 One can evaluate the quantity \R(t, x)\ = \x — y(s(x)) \ at a point near the worldline, x « y(t), solving 
(2.4). One obtains, 

, _m I v cos a + V 1 — v 2 sin 2 a \ , _ ( . . . _, 

\R(t,x)\= I f-^ 2 J \x-y{t)\+o(\x-y{t)\ 2 ), 

where v is the velocity at t, and a is the angle between v and x — y(t). Therefore, near the worldline |_R| 
represents indeed the spatial distance from the particle, a part from a never vanishing constant. 
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Moreover, away from the worldline, i.e. for x ^ y(t), we have the pointwise limit, 



hmTr(x) = rZ(x). (3.3) 

But this limit does not exist in the topology of S', due to the reemerging singularities of 
rZ along the worldline. Before taking the iS'-limit one must isolate, and subtract, these 
singular terms. 

In the present framework the form of these singular terms - and this is one more 
advantage of our regularization - is extremely simple. Since the singularities appear only 
along the wordline, the counterterm must be covariant, symmetric, traceless and supported 
along the wordline. These requirements fix it essentially uniquely, a part from an overall 
constant. We propose as renormalized energy-momentum tensor for the electromagnetic 
field generated by a point-particle, the expression, 



TZ = S'- lim 



TZ - — f (u»u u -\rf v ) 5 A {x - y{s)) ds 



= S'-\imTZ. (3.4) 



32eJ V 4 

The first counterterm, proportional to u^u v , can be viewed as an electromagnetic mass 
renormalization and is well known, see [21]. The second counterterm, proportional to rf v , 
to our knowledge has never been noticed before in the literature, probably because it can 
not be interpreted as a mass term and amounts, therefore, to a type of singularity that is 
not present in the Lorentz-equation 7 . This term is, however, necessarily present because 
the divergent part of a traceless tensor must be traceless. 

The properties of the tensor (3.4) are indeed summarized by the following 

Theorem I 

1) The limit in (3.4) exists and hence TZ £ S' ■ 

2) TZ is a Lorentz-covariant, symmetric and traceless tensor field. 

3) TZ( X ) = T em( x ) f or x ^ n the complement of the particle's wordline. 

4 ) The four-divergence of TZ amounts to, 

d,TZ = ~ J + w'uA 5\x - y{s)) ds, (3.5) 



7 The technical reason for the omission of this term stems from the fact that in the literature the 
momentum conservation equation is derived integrating the formal (divergent) energy-momentum tensor 
over a restricted class of spacelike surfaces, that are orthogonal to the four velocity u^. In this way 
one looses a divergent term proportional to - / {u^u" — i]^) S 4 (x — y(s)) ds. 
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for an arbitrary smooth worldline y^(s). 

5) If one defines the energy momentum tensor of the particle in a standard way as, 

T£ v = mj u»u v b\x - y(s))ds, (3.6) 
then the total energy-momentum tensor of the ED of a point-particle, 

= TZ + T^, (3.7) 
is conserved, d^T^ = 0, if the particle satisfies the Lorentz-Dirac equation (1.6), with 

Ft: = o. 

The crucial points to prove are the properties 1) and 4), i.e. the existence of the limit 
(3.4) and the evaluation of <9 M T e ^; we relegate these proofs to section 4. 

The limit in (3.4) holds, actually, also in a stronger sense, i.e. on a set of test functions 
that is larger then S. Indeed, since for a fixed time x° = t, for e — > the tensor T^ u (t, x) 
develops only a singularity in three-space at the point x = y(t), and due to the asymptotic 
behavior (3.2), the limit (3.4) exists also on "pseudo-test functions" of the form, 

<p(x) = 5(x° -t) <p(x), (3.8) 

where (p(x) is bounded in R 3 , and of class C°° in a neighborhood of y(t). For example, 
(p(x) can be a constant, or a characteristic function on a three-volume V. We will take 
advantage from this fact when considering momentum integrals. 

Property 2) holds by construction, and property 3) follows from the fact that the 
counterterm in (3.4) is supported entirely on the worldline. This property has its physical 
origin in the fact that the form of off the worldline is regular, and its phenomeno- 
logical consequences away from the particle, like all classical radiation phenomena, are 
experimentally very well tested. We were therefore only allowed to change the form of 
tZ on the worldline. Property 5) follows from property 4), using (1.6). 

3.1.1 Momentum integrals 

The main purpose of this subsection is to prove, using the above theorem, that the total 
four-momentum of the electromagnetic field generated by a point-particle at a generic 
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instant t is given by formula (3.16), as shown first in [8]. In the course of the proof we 
will illustrate how the new representation (3.4) can be applied in practice. 

Formula (3.4) entails indeed an operative definition for the four-momentum of the 
electromagnetic field, Py(t), contained at the instant x° = t in a volume V. We define 
this momentum applying T e ^ to the pseudo-test function, 

V (x)=5(x°-t) X v(x), (3.9) 

where Xv(%) is the characteristic function on V. By definition we have, 

p v(t) = T^tp) = hmTfV) = Jim ff(t, x) d 3 x. (3.10) 

This four-momentum is well-defined whenever y(t) ^ dV. 

If at the time t the particle is outside V, the counterterm in (3.4) does not contribute, 
and thanks to (3.3) the above definition reduces to, 

PP(t)= fr^(t,x)d\ 

coinciding with the standard expression. 

We illustrate the operative definition (3.10), computing the four-momentum Py(t) for 
a particle in uniform motion, the result being in particular useful for the derivation of 
(3.16). Due to Lorentz-invariance it is sufficient to consider a static particle in the origin. 
From (3.10) we see that we have two equivalent ways to evaluate Py(t): we can evaluate 
T^((p) on a generic pseudo-test function ip, and then set (p(x) = 5(x° — t)xv(^), or we 
can evaluate the integral f v T^°{t,x)d 3 x, and then take the limit e — > 0. For a static 
particle we are able to carry out the first, more ambitious, procedure. 

Consider first the energy density. Inserting (2.10) in (3.1) one obtains the regularized 
energy density, 

T 00 = -f—Y — r=\x\ 

2 W) (r 2 + e 2 ) 3 ' " 1 

and the 00 component of (3.4) reduces to, 

^ s '-i^(s) 2 F^-^ 3 4 
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Applying to a generic pseudo-test function ip we get, 



TZ(V) = Jim 



1 / e \ 2 r r 2 m(x° x) , d 3e 2 /■ . n , n " 

= Hm l ( ey r W£)-^,0)) d4i 
e->o 2 V47T/ 7 (r 2 + e 2 ) 3 

= 2UJ y — h (3 - n) 



where we used, 



l/e\ 2 f r 2 l3 3e 2 

d 6 x 
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2 V4W 7 (r 2 + e 2 ) 3 128 e' 

The integral in (3.11) is conditionally convergent, meaning that one has first to integrate 
over angles, and then over \x\ and x°. Inserting (3.9) we get for the energy contained in 
V the t-independent result, 

If V does not contain the origin, i.e. the particle, one has Xv(0) — and this formula 
reduces to the standard electrostatic energy. If, on the other hand, V is a sphere with 
radius p centered in the origin one has Xv(fy — 1> an d the energy in the sphere amounts 
to Py = —e 2 /8irp. In particular the total energy in whole space is zero. 

The spatial momentum P v (t) for a static particle vanishes trivially, since the Oi com- 
ponents of both terms in (3.4) are zero, even before taking e — > 0. Thanks to Lorentz- 
invariance this means that the total four-momentum of the electromagnetic field generated 
by a charged particle in uniform motion is zero, 

P"(t) = 0. (3.12) 

We consider now the total four-momentum for a particle in arbitrary motion, 

P»(t) = KmJff(t,x)d 3 x, 

where, we remember, the convergence of the integral for large is ensured by the asymp- 
totic behavior (3.2). Since for s < s the particle is in uniform motion, we have the 
additional information that, 

P M (t) = 0, for t<t, (3.13) 
14 



thanks to (3.12), where t = y°(s). 

To derive an explicit expression for P^(t) we use property 4) of the Theorem I. Define 
the vector field, 

K» = d v f v /. (3.14) 
Since the derivative is a continuous operation, (3.4) and (3.5) imply, 

S'-]imK£ = -^- [ (^ + w 2 uA5\x-y(s))ds, (3.15) 

where, again, this limit holds also on the pseudo-test functions (3.8). Integrating (3.14) 
over whole three-space we have, 

d J 7?"(f , x) d 3 x + J dif l /(t, x) d 3 x = J K£(t, x) d 3 x. 

Using the three-dimensional Gauss theorem and (3.2), the second term on the l.h.s. is 
zero. Taking then the limit e — > and using (3.15) one obtains, 

dP^(s) _e 2 (dw^ +w 2^ 
ds 6n \ ds J ' 

where we used the variable s instead of t — y°(s). Integrating this expression from — oo 
to s, and using the asymptotic relations (2.21) and (3.13), we get for the total four- 
momentum of the electromagnetic field, 

P^{s) = -—(w»(s)+ I* w 2 (X) u^(X) dX] , (3.16) 

reproducing the result of [8]. 

In appendix C we will actually prove that the more implicit expression for given 
by Rowe in [8], defines the same distribution as (3.4). 

The formula (3.16) is clearly in agreement with total four-momentum conservation, 
-P M (s) +p^(s) = constant, see (1.6). 

3.2 Covariant variational principle 

The purpose of Theorem I was the construction of a well-defined energy-momentum 
tensor. We deduce now the equations of motion and the so constructed energy-momentum 
tensor - via Noether theorem - from an action principle. 
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We propose the following regularized Lagrangian density C £ = C £ (A,y), 

C £ = -\f, v F^ - % A, -m £ J 5\x - y{s)) ds, (3.17) 

that is an element of S' if A^ is sufficiently regular. This Lagrangian differs from the 
standard Lagrangian for ED by two conceptually simple ingredients. First, the current 
j M has been replaced by its regularized counterpart see (2.11) and (2.19). Second, we 
introduced a diverging counterterm for the mass. One can indeed interpret, 

^"V. (3,8) 



m £ = m , i . 

8e UirJ ' 



as the bare mass of the regularized theory, whereas m is the physical renormalized mass. 
We define the regularized action associated to this Lagrangian in a standard way, 

I £ [A, y} = j£ £ d A x = -J Qi^F^ + j?A^ d 4 x -m £ J ds. (3.19) 

Notice that this action is gauge invariant, because the current jg is identically conserved. 
Moreover, it gives rise to the equations of motion (1.1) and (1.6) according to the following 
Theorem II 

1) The equations of motion for A 1 obtained from 1 £ are the regularized Maxwell equations, 

d,F^ = ? e . (3.20) 

These equations entail as solution the regularized LW potential A^ in (2.17). 

2) Consider the equations of motion for y M (s) derived from I £ , 

If one substitutes for A 1 the solution A^ of (3.20), one has the point-wise limit, 

ltaL MsH f _£(!£ + „v), (3.2!) 

corresponding to the Lorentz-Dirac equation. 

Property 1) is obvious. The statement 2) means that one has first to derive the 
equation of motion for y^ from I £ , then one must substitute in this equation the regularized 
LW-field potential A%, and eventually perform the limit e — > 0. This is not equivalent to 
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substituting A£ in (3.19) and then deriving the equation of motion for y^ and taking the 
limit e — > 0. This second procedure would give rise to the equation of motion ^ = 0, see 
Theorem III. The reason for the failure of this procedure is that there exists no canonical 
action giving rise to (1.6), as explained in the introduction. 

Since we are not interested in the explicit form of to prove property 2) it is sufficient 
to vary I £ w.r.t. y^, set = A^ and then take the limit e — > 0. The equation (3.21) is 
therefore equivalent to the relation, 



lim 6I e [A,y] 

£^0 



A=A £ 



= lim (-5 (m £ Jds^-J A £fl 5j? d 4 x^j 



where 5 means variation w.r.t. y^ . The proof of this relation will be given in section 5.1. 
The regularity properties of the regularized Lagrangian itself are expressed by 
Theorem III 

On the solution A^ the regularized Lagrangian converges for e — > 0, 

S'-hmC e (A e ,y) = CeS', 

where, for a generic <p G S, 

In this expression the regulator has been everywhere removed. 

The first integral in (3.23) is conditionally convergent, as (3.11). The regularized 
lagrangian admits therefore a finite limit on the solutions of (3.20). Notice in particular 
that the divergent contribution of the bare mass in (3.19) canceled out. Along the solutions 
A 1 ^, also the action I £ itself admits a finite limit. It is obtained evaluating the limit 
Lagrangian (3.23) on the pseudo-test function <p(x) = 1, 



lim/ e L4 £ ,y] = — m / ds. 



The proof of this theorem is given in section 5.2. 
As last result we state the Noether theorem. 
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Theorem IV 

The invariance of the distribution C £ under space-time translations, i. e. the identity, 

StC ° " 6A " 6 S - 5r w» - a " d ^ = °' (3 - 24) 

with 8y^ = a M , 5 A 11 = a u d u A tJi , if evaluated on the solutions of the equations of motion 
(3.20) and (1.6), implies for e — > the conservation law d^T^ v = 0, where T^ v is defined 
in (3.4), (3.6) and (3.7). More precisely, on the solution of the equations of motion one 
has the relation, 

S' - Jim 5 t C £ = a v . (3.25) 



The significance of this theorem is clear. We want only specify the following conceptual 
point. Since C £ is an element of S' the variation in (3.24) is defined in the distributional 
sense. This means that, by definition, we have, 

(5 t £ £ )(ip) = J C £ (x + a)ip{x + a) d A x — J C £ (x)ip(x) d 4 x 

= J C £ d A x + 5 J C £ <p d 4 x, (3.26) 

where 5 indicates the variation of the integral w.r.t. A, minus the variation of the integral 
w.r.t. y. This is due to the fact that the explicit dependence on x of C £ is through the 
difference x^ — y^(s) 8 . The Proof of Theorem IV will be given in section 5.3. 

3.3 The general case 

The results presented so far generalize easily to a generic set of charged point particles 
with masses m r and charges e r , and to the presence of an external field Fj^. In this case 
the total regularized electromagnetic field is given by, 

F e U = F in +J2 F (r)e' (3-27) 

r 

8 In the present case the definition (3.26) needs only to be applied to the term multiplying m e in 
(3.17) - because all other terms are regular distributions - in which case (3.26) gives, translated back to 
symbolic notation, 

S t ^ 8\x - y(s)) ds^j = -a^ (^J S\x - y(s)) ds^j - J ~ v( s )) ds > 

that vanishes identically. In practice we can then always use directly (3.24). 
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where F{^ £ is the regularized LW field produced by the r-th particle. The electromagnetic 
energy-momentum tensor, generalizing naturally (3.4), becomes, 



TZ = S' - lim 



£^0 



(3.28) 



with T^ v given by (3.1) and (3.27). There are no divergent counterterms due to the new 
interactions, because the off-diagonal terms in T^ v are well defined distributions, with 
only 1/R 2 singularities near the worldlines, and because is supposed to be regular. 
The four-divergence of this tensor is, 



8 



u r 



5 (x-y r )ds r , 



where is the unregularized LW field of the s-th particle. This relation follows directly 
from (3.5), since the off-diagonal terms in T^ v are free from overlapping divergences, and 
on them one can simply apply the Leibnitz rule to evaluate the four-divergence. The 
total energy-momentum tensor Tj^ + Tg v ', with Tg v = J2 r m r I u^u v r 8 A (x — y r (s r ))ds r , is 
then conserved if the particles satisfy the Lorentz-Dirac equations, 



dp; 

dS r 



U r 



(3.29) 



Eventually also the regularized Lagrangian density admits the natural generalization, 



3e r 2 
128 e 



)i 



5 4 (x — y r (s r )) ds r , 



where jj? is the sum of the regularized currents of the particles. This Lagrangian gives 
rise to the Lorentz-Dirac equations (3.29), and entails via Noether theorem the conserved 
energy-momentum tensor T e ^ + Tg v . Along the equations of motion C e admits again a 
finite limit for e — > in the topology of S'. 



4 Energy— momentum tensor 

This section is devoted mainly to the proof of Theorem I. The ingredients needed are the 
same as those needed in the proofs of the other theorems. For this reason we present here 
the necessary technical tools in some detail, while we will omit them in the proofs for the 
other theorems. 
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4.1 Existence of the tensor 



em 



We start with the proof of the existence of the limit (3.4), meaning that we have to isolate 
in T^ u the terms that diverge for e — > 0, if applied to test functions ip e S. Inserting 
(2.18) in (3.1) we can rewrite the regularized energy-momentum tensor as sum of three 
terms, characterized respectively by the inverse powers 1/-R 4 , 1/R 3 and l/R 2 , 

i=i 

where we count the dimensionful regulator e as a power of R, since we have R 2 = e 2 , see 
(2.16). The explicit expressions are, 

TeT = {^) 2 ^e(- R " RV + 2 ( uR > { " RU) - £2u " uU )-^ ( 42 ) 
Tg = — ^- (2{wR)R»R v + 2{uR)R^A^ - 2e 2 u^A^) - tr, (4.3) 



An) {uRf 

—) 2 t — ~~tt {-A 2 R»R V - e 2 A^A u ) - tr, (4.4) 
,4tt/ (wi?) 6 v / 

where — tr means the subtraction of the trace of the tensor, and we defined symmetrization 

by aH^ = \{a^h v + a u b> 1 ). 

We will now show that the divergent parts of these tensors, for e — > in the topology 

of S', amount to, 

T%\ div = - [£) ^lne J u^5\x-y{s))ds, (4.6) 
T e3\ d i v = 0. (4.7) 



Summing up one sees that the logarithmic divergences cancel, and that the polar diver- 
gence is canceled by the counterterm in (3.4). 

To derive (4.5)-(4.7) one must apply to a test function tp{x), and analyze the limit 
e — > 0. In doing this we encounter the technical difficulty that the kinematical variables 
y, u and w appearing in depend on x in a rather complicated way, since they are 
evaluated at the proper time s e (x). To disentangle this implicit dependence we use the 
following identity for a generic function / e C°°(R), 

f(s £ (x))= J5(s-s £ (x))f(s)ds = j '25 + [(x-y(s)) 2 -e 2 ]u»(s) (x-y(s))J(s) ds. (4.8) 
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When we apply such a function to a test function it is then also convenient to perform in 
the resulting integral the shift, 

In the following we make a systematic use of this strategy, beginning with the proof 
of (4.5). Applying (4.2) to a test function (p one obtains then, 

T?i(tfi) = A» u (tp) + B^(tp), (4.9) 

where, setting momentarily e/Air = 1, 

2bAx 2 -e 2 



Af((p) = [ds [ d A x + y e ' (-x»x v + 2{ux)u^x u) - eW - tr) (p(y(s)), 

x 2 - 

ux) 1 



B^(ip) = J ds J d 4 x 26+ \ X \ 5 (-^x" + 2{ux)u^x u) - eVu" - tr) 



Wy{s) + x)-<p(y(s)). (4.10) 

All kinematical quantities are now evaluated in the integration variable s. We have divided 
T^i((p) in two terms, adding and subtracting the term with f(y(s)). This separation is 
convenient because we will see that A^(ip) gives rise to the polar divergence in (4.6), 
while B^ u ((p) gives rise to the logarithmic divergence. 

Actually, A^ u (ip) can be evaluated exactly, it suffices to know the invariant integral, 
see appendix A, 

f d 4 x 2Mx2 ~ £2) xV = - [ d 4 x 26+ ^ 2 ~ 1} x»x» = - (5uW - rT) ■ (4.H) 
J (mi) 5 e J (mi) 5 As v ' v ' 

This gives, 

A? (if) = yj (W - \ rT) f(y(s)) ds, 

amounting to the polar divergence in (4.5). 

Consider now Bg v . From a dimensional analysis one sees that B£ u ((p) diverges loga- 
rithmically as e — > 0, i.e. as lne. To extract this divergence it is then sufficient to evaluate 
the limit, 

e^o de £ 

To compute e B^(ip) it is convenient to rescale first in (4.10) the integration variable 
x — > ex, then apply e ^, and then rescale back x — > x/e. Taking then the limit e ^ 
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the term e 2 u^u u drops out and one gets 9 , 



lim e 4- WW) = I ds I d 4 x 2S ^fJ (-x»x v + 2(ux)u^x u) - tr 
e->o de £ J J (uxf v 



•[<p(y(s)) - <p(y(s) + x)+ x a d a <p(y(s) + x)] 
3 J as 



8tt 

T 



■ J{vPw v + u v w»)(p(y(s))ds. 
This is equivalent to, 

Bf(ip) = y me j \u»w v + u v w») (p(y(s)) ds + o(l), 



where o(l) means terms that are regular for e — > 0. This gives rise to the logarithmic 
divergence in (4.5). 

The proof of (4.6) proceeds similarly, 

Tg(<p) = Jdsj d^x + | — '- (2(wx)x^x u + 2(ux)x^A u) - 2e 2 u^A v) - tr) <p(y(s)+x), 

where here A M = (ux) — (wx) u^. This is again logarithmically divergent and one 

obtains, operating as above 10 , 

d ft 2r) ( r r 2 \ 

lim e — Tg{ip) = J ds J d A x ^ ' (2{wx)x^x u + 2{ux)x^A v) - tr) x a d aV {y{s) + x) 

= -y J(u»w v + u u w») V (y{s))ds. 

This proves (4.6). 

The tensor Tf£ admits trivially a finite limit for e — > 0, since only 1/|-R| 2 powers are 

present. This concludes the proof of property 1) of the theorem. 

9 The resulting expression can be evaluated explicitly noting the integral, with y = y(s) 1 

f d 4 x H£±£J2 x p x v [<p(y) - <p(y + x )+ x a d a <p(y + x)] = - fdSl m a dMy) = 

J yllX) J yUTflj 

y {u p d v ip{y) + u v d^{y) + rTWdMy) - &uVu a d aV {v)) , 

where J dCl denotes the integral over the three dimensional solid angle, and m p = (1, jj^). 
10 This time one needs the integral, with y = y(s), 

f . 25 + (x 2 ) ON f m! 1 m v m p , , 

J {ux) 1 ' J (um) a 

^ {rf v u p + rj vp u p + rj pp u v - &u^u v u p ) <p(y), 
where J dfl is the integral over angles and m p — (1, ^). 
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4.2 Conservation of the energy— momentum tensor 

The conservation of the total energy-momentum tensor is reduced to the proof of property 
4), i.e. to the evaluation of d^TZ, see (3.5). Since the derivative is a continuous operation 
in S', once we have established the existence of the limit (3.4), the operator can be 
interchanged with the limit, and we get, 

d,TZ = S'- lim (d»Tr -^-J( w "-\ dV ) ~ v( s )) d >) ■ 
Since also the regularized LW-field satisfies the Bianchi identity, 

one deduces the standard result, 

d,Tr = (4.12) 

and hence, 

d,TZ = S>- lim [~F^j e , -JL-J(tf>-1 5\x - y{s)) ds^j . (4.13) 

The calculation is therefore reduced to the evaluation of F^j £fJl for e — > 0, in the topology 
of S'. 

A somewhat delicate point in the evaluation of this limit is the following. The formula 
for the current (2.19) carries a factor of e 2 in front, due to the fact that in the limit e — > 
the current (x) goes to zero for x^ in the complement of the wordline. Actually, it can 
be seen that when evaluating the limit (2.14), there is only one term in (2.19) that has a 
non- vanishing limit, more precisely, 

S' -lim (e 2 — 3U ^ ) = 

™ { 4tt (uR) 5 J 3 ' 

while all other terms converge to zero. These "vanishing" additional terms are however 
needed in the regularized current to ensure its conservation, and they contribute in the 
product F^j^ for e — >■ 0, since F^ has polar terms. Multiplying out (2.18) and (2.19) 
we obtain indeed, 



2 

2 



2^ 2 (bR) \ (wR) - 1 

(uR)s (uRy) (uR)* 



» "f:," , f^(M-') g 



uR) 7 V (uR) r 
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+ o(e), (4.14) 



where = dw^/ds. o(e) denotes terms that vanish for e — > in S'. In the present case 
these are terms of the form e 2 /R 4: , that converge to zero by power counting. We have now 
to perform the limit of this expression for e — > in S'. This analysis can be performed 
using the techniques developed in section 4.1, with the difference that now one has to 
keep also terms that are finite for e — > 0. The computation is deferred to appendix B 
where we show that, 



4n) 



8vr (dw v 2 \ 7T 2 / 1 ot 



3 V ds 2e\ 4 



6\x-y(s))ds + o(e), (4.15) 



where o(e) denotes again terms that vanish for e — > in S' . 

Inserting this in (4.13) one sees that the polar contributions cancel, as they must, and 
one obtains (3.5). 

4.3 Comparison with Rowe's energy— momentum tensor 

The electromagnetic energy-momentum tensor introduced by Rowe [8] is written as the 
sum of three terms elements of S', 

©£ = E©r, ( 4 - 16 ) 

i=l 



where, 



0f = d a Kr + ^ [u^5 A (x-y(s))ds, (4.17) 

0f = d a Kr-^-JuV5\x-y( S ))ds, (4.18) 
/ e \ 2 A 2 

e f = -It-) tt^ w - ( 4 - 19 ) 

6 \AttJ (uR) 6 v 1 

The tensors K^ v e 5' are antisymmetric in a and /x, and are given by, 

—) -^—(R^A a -R a A^)R v . (4.21) 

An) (uRf K ' K ' 



K. 



In these expressions no regularization is needed. Notice, indeed, that K^" contains only 
1/i? 2 singularities and that K"^ contains only derivatives of 1/R 2 singularities, so that 
both tensors are elements of S'. This implies that also 0^ and 02 u are elements of 
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S', because the derivative of a distribution is again a distribution. The terms with the 
5-functions in (4.17), (4.18) are added to ensure the symmetry of each Of separately - 
a property that is rather hidden in Rowe's formulation, and needs to be proven a poste- 
riori. A more serious practical drawback of the expressions above is that the derivatives 
appearing are distributional derivatives, and they can not be evaluated applying simply 
the Leibnitz rule, due to the singularities present on the worldline. 

Away from the worldline, i.e. for x M 7^ y M (A), one can evaluate the derivatives above 
using the Leibnitz rule, and it is then straightforward to check that in the complement of 
the worldline coincides with the sum of (4.2)-(4.4), with e — 0, and hence with our 
TZ in (3.4). 

The comparison of 0^ m with T e ^ as elements of S' can be performed applying to 
the tensors K^ u our iS'-regularization, K^ u — > K^ u , where the regularized tensors are 
given formally again by (4.20), (4.21), but with the replacement s(x) — > s £ (x). On the 
regularized tensors one can now use the Leibnitz rule to evaluate d a K"j lu , and perform 
eventually the <S'-limit for e — > 0. Following this procedure in appendix C we prove indeed 
that, 

C\^V rpfiv 

em em' 

as distributions. Here we note only that a comparison of (4.19) with (4.4) yields immedi- 
ately, since there are only 1/R 2 singularities present, 

S'-hmT^ = Of. (4.22) 



5 Proof of variational results 

5.1 Derivation of the Dirac— Lorentz equation 

In this subsection we proof property 2) of Theorem II, that is equivalent to the relation 
(3.22). Instead of computing the variation of I £ under a generic variation of we 
compute the variation of the Lagrange density C £ . The limit of SI £ for e — > can then 
be obtained integrating 5C £ between two hypersurfaces, and sending e to zero. Another 
reason for proceeding in this way is that the explicit form of SC £ will also be used in the 
proof of Theorem IV. 
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In this subsection with "<5" we will always mean a generic smooth variation w.r.t. y^. 
Starting from (3.17) we have then, 

5C £ = -A £ ,5j? - m £ 5 (/ 5\x - y(s)) ds^j , (5.1) 

where, according to what we need in (3.22), we have replaced - after variation - — > A £ll . 
While the variation of the second term is trivial, to vary the first term it is more convenient 
to write, 

-A^5j» = -5 (A £ ^) + 5A £ ,j», (5.2) 

since SA £fJi is simpler then 8j%. The variation of the first term in this expression can be 
obtained from the ^'-expansion, 

j?A»= j (^uW-47ruW + o(e)y 4 (x-y( S ))ds, (5.3) 

upon contracting /x and v. This expansion can be derived using the method illustrated in 
appendix B, see (7.7). 

The evaluation of the second term in (5.2) is more cumbersome, and the calculation 
is deferred to appendix D. The result is, 

SA ~* = (sH/(-l>«-f + 

+<9 M | (^7r(w5y)u^ -^(5y^-(u5y)u^5 4 (x-y(s))ds +o(e).(5.4) 

Using this in (5.2) one obtains eventually, 

SC £ = J (mm* - ^ + w 2 u^j ^ % 5\x - y{s)) ds + (5.5) 

+d, J [^{w 8y) u»+(m+ (5y» - (u 5y)u^ 5\x - y{s)) ds + o(e), 

which is our final result, holding for generic variations of y. The first line is finite as e — > 0, 
and the second line, which corresponds to a four-divergence, contains a polar term in e. 

To derive the equations of motion we compute 5I £ integrating 5C £ between two hy- 
persurfaces, and impose that on them 5y^(s) vanishes. The four-divergence drops then 
out - together with the polar term - and the first line (5.5) gives (3.22). This concludes 
the proof of Theorem II. 
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5.2 Proof of Theorem III 



The evaluation of the regularized Lagrangian on the solution = A^ requires to evaluate, 
using (2.18) and (7.6), 




This gives, using (5.3), 

C C {A„V) = -jWi'" -A c ^-m c js\x- y (s))ds 

Applying to a test function one gets, 

Km£ £ (A £ ,y)( V ) = jjm j " ds j ™ +{ fZ ^ [<p(x + y(s)) - <p(y(s))] d*x 



where we used, 




The integral above converges now for e — > 0, and we can set e = obtaining the condi- 
tionally convergent integral, 

1 / e \ 2 r f 25 (x 2 ) t 
km£ £ (A £ ,y)(ip) = -[—) J dsj -^L[ ip { x + y{ s ))- V {y{s))}d 4 x-mJ V {y{s))ds. 

Shifting x — > x — y(s) and integrating out the 5-function one obtains (3.23). 

5.3 Proof of the Noether Theorem IV 

The translation invariance of the regularized Lagrangian is expressed by the identity 
(3.24), 5 t C £ = 0, where C £ is given in (3.17), and Sy* = < 5A» = a v d v A». 
After standard steps one obtains, 

5 t C £ = (d,F^ - f £ )5A u - d,(5A u Fn + A^ £ + S y (m £ J 5\x - y(s)) ds^j - a^C £ , 
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where 5 y means variation w.r.t. y p . We evaluate this expression now on the solutions of 
the regularized Maxwell equation, A p = A p , and on those of the Lorentz-Dirac equation. 
The first term goes then to zero, and for the combination, 



Ae^yf e + $y (m £ J 5 4 (x - y(s)) ds^j , 



we can use the general expansion (5.1), (5.5), with <5y M = a M , with opposite sign. Thanks 
to the Lorentz-Dirac equation, the expression (5.5) reduces to a four-divergence and hence 
5 t C e becomes a four-divergence, too. Using again (5.3) to evaluate A £U j" appearing in 
<9 M £ £ , after simple algebra one obtains, 

6 t £ e = a v d li fr + o(e) = 0, (5.6) 

where, 

Tf = F p /d v A ep + - r]^F e pa F epa + m J u p u u 5\x - y(s)) ds 

+ (^) 2 / u»u v - 4ttu p w u + ^ rTj 5\x - y{s)) + o(e). (5.7) 

Even if it is not obvious from the derivation - see however the footnote below - it can be 
seen hat the limit, 

S'-\imfr = f^, (5.8) 

exists. (5.6) implies then that the energy-momentum tensor T pv is conserved. It is 
immediately seen, however, that this tensor does not coincide with the one constructed 
in Theorem I, see (3.7), nor is it symmetric in \x and v. On the other hand this was to be 
expected, since the Noether theorem leads in general to the canonical energy-momentum 
tensor, not to the symmetric one. 

From the general theory we know however that, if the Lagrangian density is not only 
translation invariant but also Lorentz invariant, there exists always a symmetrization 
procedure for the canonical tensor. In the case of classical Electrodynamics this procedure 
requires to add the divergenceless term —d p (F pp A u ). In the case at hand this suggests 
then to add the divergenceless term 11 , see (5.3), 

11 The product F^A v e contains terms that behave as l/R 3 , if e = 0. This means that a priori its 
(S'-limit for e — > contains logarithmic divergences, ~ Ins. However, by inspection i.e. applying the 
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= -F p /d p A v e + J (^-^uV + 4ttu^w u + o(e) ] j5 4 (x-y(s))ds. 

Adding this to (5.7) one sees that eventually one gets a symmetric tensor, being indeed, 

ff - d p {F^A v e ) = T» v + o(e), (5.9) 
where T^ v is defined in (3.7). This concludes the proof of the Noether theorem. 
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6 Appendix A: Invariant integrals 



In the proofs throughout the paper one encounters the invariant integrals, 

_ r 4 25 + (x 2 - 1) 

■*■ JJ, J CI X / \ AT X 7, 



I fxv i d x , \ at XuXl 





{ux) N 






_(x 2 - 


1) 




(ux) N 
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where S + (x 2 — 1) = H(x°)S(x 2 — 1). These integrals can be evaluated considering ■u' 1 as 
an unconstrained variable, not satisfying u 2 = 1, and writing, for example, 
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This reduces these integrals to the calculation of the "generating function" , 

25 Ax 2 -i) vr 3 / 2 r(f-i; 
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product to a test-function and taking e — > one checks that these divergences cancel out, and the limit 
exists. This means that also, 

5'-lim8 f (TO 

e— >0 

exists. This fact, together with (5.9), provides an indirect proof for the existence of the limit (5.8). 
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After taking the derivatives w.r.t. one sets again u 2 = 1. For example, 
and hence, 



(mx) 3 (w 2 ) 3 / 2 ' 
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Setting -u 2 = 1 this gives (4.11). 



7 Appendix B: proof of (4.15) 

Before taking the limit of e — > it is convenient to rewrite (4.14) in the form, 
2 // 2w 2 
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We have used, 
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and we have omitted the term, 
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,{wR)w l 



{uRf ' 

because it goes to zero for e — > in 5', for kinematical reasons. More precisely, after 
applying it to a test function and integrating over d^x, for invariance reasons the factor 
R^ in the numerator gets replaced by w M , and (uw) = 0. This will become clear below. 
Regarding the last term in (7.1) it is convenient to note the identity, 
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+ o(e), (7.2) 



(uR) 7 (uR) e 

where the o(e) terms are in the same sense as above. Inserting this in (7.1) one sees that 
the (W?)-term cancels getting, 
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The behaviour of the three terms present, for e — > in S', are as follows, 

^ R v = y / ^«^ 4 (x - y(a)) da + (e), (7.4) 
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For the proofs of these relations, that require to apply the l.h.s. to a test function and 
to analyze the limit e — > 0, one can use the technicalities developed in section 4.1. We 
report explicitly the proof of (7.5), 
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where we used the invariant integrals of appendix A. 
Inserting (7.4)-(7.6) in (7.3), one obtains (4.15). 

8 Appendix C: Rowe's energy— momentum tensor 

Before comparing of (3.4) with 6^ in (4.16), we cast the former in a slightly different 
form writing, see (4.2)-(4.4), 
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where, 



31 



The advantage of this form of presenting T e ^ is that the tensors converge separately 
for e — > in the topology of S', as does . We have added and subtracted the same 
term from and , to eliminate the logarithmic divergence from both. This can be 
seen from (4.5)-(4.7), noting that, 
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We evaluate now the four-divergencies appearing in (4.17), (4.18), using the regular- 
ized versions K°f ,v for (4.20), (4.21), obtained replacing s(x) — > s £ (x). Then one can 
apply Leibnitz' rule and, taking into account that now R 2 = e 2 ^ 0, after a straightfor- 
ward computation one obtains, 
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In d a K^i U we have put in evidence the polar part, proportional to u^u v — \r\^ v , and in 
d a K°T we have set = dw^/ds. 

We can now consider the <S'-limit for e — > of these relations, applying them to a 
test-function and sending e — > 0. The calculations, that are straightforward but a bit 
lengthy, can be performed with the techniques illustrated in appendix B, see (7.7). Since 
the limits of d a K^ v for e — > exist by construction, the r.h.s. of the above relations must 
also admit finite limits. It can indeed be seen that the square bracket in (8.3) exhibits 
a divergent contribution that cancels exactly the polar term in (8.2), while the square 
bracket in (8.4) admits a finite limit by power counting. The results are, 

8 a K^ v = (s'~ limDg) + — I (\u^w v + u u w») - — u^dA 5\x - y(s)) ds, 

\ e^O J 7T J \6 16 / 

d a KT v = (s'-hmD^-^-Ju»w»6\x-y(s))ds. 
Inserting these expressions in (4.17), (4.18) one gets the symmetric and traceless tensors, 
0f = (s' - Jim^r) + ^ J (uV + u u w^)5\x - y{s)) ds, (8.5) 
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0f = (S'~ lim D%) - —I {u»w v + u u w»)5\x - y{s)) ds. (8.6) 
These formulae, together with (4.22) and (8.1), imply then 0^ = T e ^, q.e.d. 



9 Appendix D: Proof of (5.4) 



To evaluation of 5A £ll j£ requires first to evaluate SA £/I under a generic smooth variation 
of y* 1 . This evaluation is complicated by the fact that A^ - see formula (2.7) - depends on 
y^ explicitly, but also implicitly through the (regularized) retarded proper time function 
s £ (x), defined in (2.8). To determine 5A £fM it is convenient to use the following shortcut. 

We compute 5A^ applied to a test function ip e S, inserting a 5-function as in (4.8), 
and then shifting x — > x + y(s), 
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where in the square bracket all kinematical quantities are evaluated at s = s e (x). This 
gives, 
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where, 
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The appearance of the combination £ M is due to the invariance of A^ under reparametriza- 
tion of the worldline, i.e. 5A^ = if 5y^ = 5\u^. 
Multiplying out (9.1) and (2.19) one obtains then, 
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where the terms o(e) correspond to expressions of the type e 2 /R 4 , that converge to in 
the topology of S', by power counting. To evaluate the limit of this expression under 
e — > in S', one proceeds precisely as in (7.7). In this way the evaluation of the limit of 
all terms above is reduced to the determination of the invariant integrals of appendix A. 

Lets consider e.g. in the expression above the most singular term, corresponding to 
e 2 (R£) I (uR) 8 , that by power counting would give rise to a double pole 1/e 2 in that, 
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We compute the invariant integrals according to appendix A, 
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Since (w£) = we see that the double pole cancels, and the remaining two terms give a 
1/e contribution and a finite part. After integration by parts one obtains, 



_ 2 m 



i.e. in S', 



All other terms in (9.2) by power counting give simple pole contributions and finite parts, 
too, and a straightforward but a bit lengthy calculation gives (5.4). 
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